Multiband Transport in Bilayer Graphene at High Carrier Densities 
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We report a multiband transport study of bilayer graphene at high carrier densities. Employing 
a poly (ethylene) oxide- CsC104 solid polymer electrolyte gate we demonstrate the filling of the high 
energy subbands in bilayer graphene samples at carrier densities \n\ > 2.4 x 10^^ cm~^. We observe 
a sudden increase of resistance and the onset of a second family of Shubnikov de Haas (SdH) oscilla- 
tions as these high energy subbands are populated. From simultaneous Hall and magnetoresistance 
measurements together with SdH oscillations in the multiband conduction regime, we deduce the 
carrier densities and mobilities for the higher energy bands separately and find the mobilities to be 
at least a factor of two higher than those in the low energy bands. 

PACS numbers: 73.63.b, 73.22.f, 73.23.b 



Multiband transport is common for many complex 
metals where different types of carriers on different pieces 
of the Fermi Surface (FS) carry electrical currents. Con- 
duction in this regime is controlled by the properties of 
the individual subbands, each of which can have distinct 
mobilities, band masses, and carrier densities. Other 
changes to the single-band conduction model include 
inter-band scattering processes and mutual electrostatic 
screening of carriers in different subbands, which alters 
the effective strength of the Coulomb potential and hence 
adjusts the strength of electron-electron and electron- 
charged impurity interactions. 

To understand electronic conduction in this regime, 
it is desirable to study the properties of the individ- 
ual bands separately and compare these to the prop- 
erties in the multiband regime. This was achieved in 
2-dimensional electron gases (2DEGs) formed in GaAs 
quantum wells [1], where the subbands can be contin- 
uously populated and depopulated by inducing parallel 
magnetic fields. In these 2DEGs, an increased overall 
scattering rate due to interband scattering was observed 
upon the single- to multiband transition, [I, along 
with changes in the effective Coulomb potential which 
led to the observation of new filling factors in the frac- 
tional quantum Hall effect M. 

Bilayer graphene (BLG)[5|49|, with its multiband struc- 
ture and strong electrostatic tunability, offers a unique 
model system to investigate multiple band transport phe- 
nomena. BLG's four- atom unit cell yields a band struc- 
ture described by a pair of low energy subbands (LESs) 
touching at the charge neutrality point (CNP) and a 
pair of high energy subbands (HESs) whose onset is 
- ±0.4 eV away from the CNP (Fig.l(a)). Specifica%, 
the tight binding model yields the energy dispersion [7|: 
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FIG. 1: (a) The tight-binding band structure of bilayer 
graphene for interlayer asymmetries A = eV (gray) and 
A = 0.6 eV (black), (b) Schematic view of the double gated 
device, consisting of the Si02/Si back gate and the electrolytic 
top gate. Debye layers of Cs^ or CIO4 ions are formed 
d ^ 1 nm above the bilayer and the gate electrode, respec- 
tively, (c) Longitudinal resistivity and Hall resistance of the 
bilayer graphene device at T = 2 K as a function of back gate 
voltage Vbg for 3 different fixed electrolyte gate voltages Veg = 
-1.7, -0.4, and 1 V from left to right, corresponding to pre- 
doping levels of uh — (-2.9, 0, 2.9) x 10^^ cm~^. Inset shows 
an optical microscope image of a typical Hall Bar device (the 
scale bar corresponds to 5 /xm). 



is the wave vector measured from the Brillouin zone cor- 
ner, vf ^10^ m/s is the Fermi velocity in single layer 
graphene, 71 ^ 0.4 eV is the interlayer binding energy, 
and A is the interlayer potential asymmetry. Interest- 
ingly, since a perpendicular electric field E across the 
sample gives rise to an interlayer potential difference A, 
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FIG. 2: Landau fan diagram of the differential longitudinal resistivity dpxx/duH for 3 different density ranges at T = 2 K 
as a function of the Hall density uh and the magnetic field, (center) The SdH oscillations in the LES converge at the CNP 
and flatten out at higher uh due to decreasing LL separation, (left and right) For |n//| > 2.6 x 10^^ cm~^ additional SdH 
oscillations appear, originating at the resistivity spikes (red regions) that mark the onset of the HES. 



it opens up a gap in the spectrum of the LES H, [lo|, [ill 
and is furthermore predicted to adjust the onset en- 
ergy of the HES. Whereas the LESs have been widely 
studied, the HESs, with their expected onset density of 
n > 2.4 X 10^^ cm~^ have thus far not been ac- 
cessed in transport experiments. This can mainly be 
attributed to the carrier density limitations set by the 
dielectric breakdown of the conventional Si02/Si back 
gates, which do not permit the tuning of carrier densities 
above n ^ 0.7 x 10^^ cm-^ (e^ ^ 0.2 eV). 

In this letter, we report multiband transport in bilayer 
graphene. Using an electrolytic gate, we were able to 
populate the HES of bilayer graphene, allowing for both 
the LES and HES to be occupied simultaneously. The 
onset of these subbands is marked by an abrupt increase 
of the sample resistivity, most likely due to the opening 
of an interband scattering channel, along with the ap- 
pearance of a new family of Shubnikov-de Haas (SdH) 
oscillations associated with the HES. A detailed analy- 
sis of the magneto- and Hall resistivities in combination 
with the HES SdH oscillations in this regime enables us 
to estimate the carrier mobilities in each subband sepa- 
rately, where we observe a two-fold enhanced mobility of 
the HES carriers as compared to the LES carriers at the 
same band densities. 

Bilayer graphene devices were fabricated by mechani- 
cal exfoliation of Kish graphite onto 300 nm thick Si02 
substrates, which are backed by degenerately doped Si 
to form a back gate. The samples were etched into a 
Hall bar shape with a typical channel size of ~ 5 jiia 
and then contacted with Cr/Au (0.5/30 nm) electrodes 
through beam lithography (Fig. 1(c) inset). In order 
to access the HES we utilized a recently developed solid 
polymer electrolyte gating technique [Ill4l4j. which was 
recently shown to induce carrier densities beyond values 
of n > 10^^ cm~^ [14] in single layer graphene. The 
working principle of the solid polymer electrolyte gate is 



shown in Fig. 1(b). Cs+ and CIO4 ions are mobile in the 
solid matrix formed by the polymer poly (ethylene) oxide 
(PEO). Upon applying a gate voltage Veg to the elec- 
trolyte gate electrode, the ions form a thin Debye layer 
a distance d ^ 1 nm away from the graphene surface. 
The proximity of these layers to the graphene surface re- 
sults in huge capacitances per unit area Cg^, enabling 
extremely high carrier densities in the samples. While 
CSCIO4 has almost the same properties as the typically 
used LiC104 salt, we find a reduced sample degradation 
upon application of the electrolyte on top of the sample, 
resulting in considerably higher sample mobilities. 

One major drawback of the electrolyte gate for low 
temperature studies is that it cannot be tuned below 
T < 250 K, where the ions start to freeze out in the 
polymer and become immobile (though leaving the De- 
bye layers on the bilayer surface intact) [l2|, [ij]. A de- 
tailed study of the density dependent transport prop- 
erties at low temperatures can therefore be quite chal- 
lenging. In order to overcome this issue, we employ the 
electrolyte gate just to coarsely tune the density to high 
values (|n| <10^^ cm-^) at T = 300 K, followed by an 
immediate cooldown to T = 2 K. We then use the stan- 
dard Si02/Si back gate to map out the detailed density 
dependence of the longitudinal sheet resistivity pxx and 
the Hall resistance Rh-, from which we extract the total 
carrier density of the sample uh = B/eRn, with B the 
magnetic field and e the electron charge. Here we find 
the back gate capacitance to be C^g = 141 aF//im^, al- 
most unaltered by the presence of the Debye layers on 
top of the sample. 

In this experiment, we have measured pxx and Rh of 
more than 3 BLG devices as a function of the back gate 
voltage Vbg at various fixed Veg corresponding to the 
wide density range of uh ±8 x 10^^ cm~^. Fig. 1(c) 
shows Pxx and Rh for a representative device for 3 se- 
lected cool-downs at Veg = -1.7, -0.4, 1 V from left to 
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right, corresponding to a pre-doping level of uh =(- 
2.9, 0, 2.9) X 10^^ For low doping levels {Veg =- 

0.4 V, Fig. 1(c) middle) we observe the expected Dirac 
Peak in p^^ and the ambipolar transition of Rh as V^g 
sweeps through the CNP. Away from the CNF, pxx and 
Rh decrease as \tih\ increases, as was observed before 
in BLG samples [9]. For the strongly pre-doped gate 
sweeps however (Fig. 1(c) left and right), we observe a 
rather unexpected non-monotonic feature in the sam- 
ple resistivity. Instead of a monotonic decrease of pxx 
with increasing {uhI^ it exhibits an abrupt increase by 
~ 10% symmetrically at both electron and hole sides at 
n* ~ In// 1 = 2.6 x 10^^ cm~^, a carrier density which 
is consistent with theoretical expectations for the onset 
density of the HES 0, HI- A similarly increasing resis- 
tivity at the opening of a new subband was previously 
observed in 2D electron systems formed in wide GaAs 
quantum wells [l|-[3|, where it was attributed to a de- 
creased overall scattering time r due to the opening of 
an additional inter-band scattering channel as the new 
subbands are populated. Such an inter-band scattering 
mechanism between the LES and the HES is also ex- 
pected to give rise to a resistivity increase upon filling 
of the HES in BLG samples, making it a likely candi- 
date for the origin of the observed resistivity increase. 
However, considering the strong differences between the 
2D electron gases in GaAs quantum wells and in BLG, 
including vastly different densities of states, mobilities, 
and electron energies, more theoretical work needs to be 
done to conclusively determine the origin of this resistiv- 
ity increase. 

The electronic structure of the LES and HES can be 
further investigated by studying the effect of the mag- 
netic field B on the longitudinal resistivity pxx{B) in the 
various density ranges. Fig. 2 shows the Landau fan di- 
agram of the differential sheet resistivity dpxx/<^'riH as a 
function of B and uh- Close to the CNP (Fig. 2 center) 
the SdH oscillations in the two LES are quite pronounced, 
but with increasing density their amplitude quickly de- 
cays as the energy separation of the Landau Levels (LL) 
decreases. Above the onset of the HES (Fig. 2 left and 
right), marked by the "spikes" of increased resistivity 
(here the red regions) however, we observe another set of 
SdH oscillations which form LL fans converging into the 
onset point of the HES. 

In order to analyze the SdH oscillations, we now plot 
the pxx{B) traces for various fixed uh as a function of the 
inverse magnetic field B~^ . Fig. 3(a) displays three ex- 
emplary traces at different uh above the onset density of 
the HES. All traces show periodic oscillations in B~^ al- 
lowing us to obtain the SdH density, nsdH = x ^(^~^)' 
assuming that each LL is both spin and valley degen- 
erate. Whereas for all \nH\ < n* we find that the ob- 
tained nsdH ^ TiH^ indicating that the SdH oscillations 
are solely from a single band (i.e., the LES), for \nH\ > n* 
the obtained nsdH values are much smaller than the si- 
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FIG. 3: (a) Exemplary traces of the longitudinal resistivity as 
a function of inverse magnetic field at fixed values of uh be- 
yond the onset of the HES. (b) Carrier densities inferred from 
the SdH oscillations vs. the overall Hall densities uh^ from 4 
cooldowns at different set electrolyte gate voltages Veg — -2Y 
(red), -1.7 V (orange), 1 V (purple), 1.4 V (blue), (bottom) 
\nHEs\ vs. riH^ fitted with theoretical expectations for the 
HES. (top) \nLEs\ vs. uh^ fitted with theoretical expecta- 
tions for the LES. Line traces correspond to theoretical fits 
for different values of A = 0.31 eV (red), 0.17 eV (orange), 
0.13 eV (purple), 0.26 eV (blue). 

multaneously measured uh values. This behavior can be 
well explained by assuming that these SdH oscillations 
refiect only the small fraction of charge carriers lying in 
the HES. For \nH\ > ^* we hence are able to extract 
the occupation densities of the LES (tiles) and HES 
(riHEs) from riLES = riH - usdH and uhes = nsdH- 
Fig. 3(b) shows the \nLEs\ and \nHEs\ in this regime 
as a function of the total carrier density \nH\' For each 
fixed Veg^ the obtained \nLEs\ and \nHEs\ increase as 
\nH\ increases (adjusted by V^g)^ for both electrons and 
holes. Interestingly, we notice that the \nLEs{'^H)\ are 
slightly larger for larger | Veg \ while the trend is opposite 
for the HES, i.e. \nHEs{'^H)\ are smaller for larger \Veg\^ 
even though their uh values are in similar ranges. These 
general trends can be explained by an increase of the 
interlayer potential difference A for increased values of 
\Veg\^ which are predicted by the tight-binding model in 
Eq.l to result in an increase of the onset density (energy) 
of the HES. 

While a precise quantitative determination of the ex- 
pected shift in the onset density of the HES as a func- 
tion of Veg and V^g requires a self-consistent calculation 
of /S^iVeg-, Vhg) and would go beyond the scope of this pa- 
per, we can still qualitatively test the above prediction. 
This is possible since A is mostly controlled by Veg , which 
has a much stronger coupling to the BLG sample than 
the Vhg^ thus allowing us to approximately treat A as a 
constant for fixed Veg • Since the experimental traces dis- 
played in Fig. 3(b) correspond to different values of Veg 
but the same ranges of V^g^ A is different for each trace 
and can be extracted from the theoretical fits from Eq. 1, 
with A as the only fitting parameter. Indeed for all 4 
traces we find good agreement with the theoretical fits; 
we clearly observe an enhanced onset density (energy) 
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for the traces with larger set potential differences across 
the sample, which is in good qualitative agreement with 
theoretical predictions. 

We now turn our attention to the transport properties 
of BLG in the limit of uh > tl* . The filling of these 
sub-bands creates a parallel transport channel in addi- 
tion to the one in the LES, thus defining the transport 
properties in this regime by two types of carriers with 
distinct mobilities /ii,2, effective masses 2 and sub- 
band densities ni,2 (here the index corresponds to the 
LES(l) and HES(2)) [13, [Tel. In sharp contrast to a sin- 
gle band Drude model, where pxx{B) does not depend on 
the B field, in a two-carrier Drude theory it is expected 
to become strongly modified, resulting in a pronounced 
B field dependence [l7j : 



.4 B(T) , 



Pxx{B) 



ni/ii + n2/i2 + (^iMiM2 + n2/i2Mi)^^ 
e{{nipi + n2/i2)^ + lAlAi^i + ^2)^^^) ' 



(2) 



Fig. 4 (a) shows magnetoresistance traces for different 
fixed Hall densities uh- Close to the CNP, where only 
the LES are populated (Fig. 4(a) black trace), the pxx{B) 
traces are nearly fiat as expected from the one-fiuid 
Drude theory. When the density is increased and the 
HES starts to fill up, however, we observe a smooth tran- 
sition to an approximately parabolic B field dependence, 
resulting in a strong increase of pxx of up to 25% from 
T to 8 T. Using the previously extracted carrier den- 
sities in the two bands ni^2 we can now fit the pxx{B) 
traces with the two-carrier Drude model in Eq.2, with 
the mobilities of the two subbands pi^2 as the only fit- 
ting parameters. As shown in Fig. 4(b) the experimental 
finding are in excellent agreement with the theory, allow- 
ing us to deduce the values of pi^2 with good accuracy. 
Moreover, the ability to extract the mobilities of the HES 
allows us now to characterize the HES in more detail. 

In Fig. 4(c) we plot the extracted mobilities of the HES 
P2 against the carrier density in the HES 712 and compare 
it to the mobilities pi of the LES at a similar range of 
subband densities in the LES ni. We find that the mo- 
bilities in the HES are at least a factor of two higher than 
those in the LES. Considering that the effective carrier 
masses are similar for the LES and the HES, this fea- 
ture of the HES may be due to the enhanced screening 
of charged impurity scatterers at higher carrier densities, 
effectively reducing the scattering rate of the HES carri- 
ers on these scatterers. A more detailed theoretical study 
is required, however, to undertake a quantitative analysis 
of this problem. 

In conclusion, using a polymer electrolyte gate we have 
achieved two-band conduction in bilayer graphene. We 
have found that the filling of these bands above a Hall 
density of {nnl > 2.4 x 10^^ cm~^ is marked by an in- 
crease of the sample resistivity by ~ 10% along with the 
onset of SdH oscillations. From simultaneous Hall and 
magnetoresistivity measurements, as well as the analy- 
sis of the SdH oscillations in the two carrier conduction 
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FIG. 4: (a) Longitudinal resistivity pxx{B) as a function 
of magnetic field for uh = (0.84, 2.88, 3.00, 3.31, 3.83, 
4.21, 6.96) X 10^^ cm"^ from top to bottom. The pxx{B) 
traces undergo a smooth transition from a nearly B indepen- 
dent behavior when the HES is empty (below uh < n* ^ 
2.6 X 10^^ cm~^), to a strong, non-trivial B dependence when 
the HES is occupied, (b) An exemplary pxx{B) trace at 
riH =6.96x10^^ cm"^ and nsdH =1.42x10^^ cm"^ with ac- 
companying fit (dashed line) from Eq.2, using the mobilities 
in the LES pi = 541 Vs/cm^ and the HES p2 = 2428 Vs/cm^ 
as fitting parameters, (c) The mobilities /ii,2(n) as extracted 
from pxx{B) traces at various fixed uh as a function of the 
density in the individual subbands. 



regime, we have characterized the distinct carrier densi- 
ties and mobilities of the individual subbands, where we 
have found a strongly enhanced carrier mobility in the 
HES of bilayer graphene. 
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